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Abstract Mountain lee waves have been previously observed in data from the Moderate Resolution 
Imaging Spectroradiometer (MODIS) "water vapor" 6.7 pm channel which has a typical peak sensitivity at 
550 hPa in the free troposphere. This paper reports the first observation of mountain waves generated by the 
Appalachian Mountains in the MODIS total column water vapor (CWV) product derived from near-infrared 
(NIR) (0.94 pm) measurements, which indicate perturbations very close to the surface. The CWV waves are 
usually observed during spring and late fall or some summer days with low to moderate CWV (below -2 cm). 
The observed lee waves display wavelengths from 3-4 to 15 km with an amplitude of variation often comparable 
to -50-70% of the total CWV. Since the bulk of atmospheric water vapor is confined to the boundary layer, this 
indicates that the impact of these waves extends deep into the boundary layer, and these may be the lowest level 
signatures of mountain lee waves presently detected by remote sensing over the land. 


1. Introduction 

Water vapor is a potent greenhouse gas the bulk of which is confined to the boundary layer. It is a major 
component of the Earth energy and water cycles and a major parameter in weather forecasting and climate 
modeling. Water vapor responds to surface temperature changes through radiative effects and climate feed- 
backs [Trenberth et al., 2005; Wagner et al., 2006; Solomon et at., 2010]. Accurate knowledge of atmospheric 
water vapor is crucial to correct the time delay and phase distortions for the repeat-pass Interferometric Synthetic 
Aperture Radar applications [e.g., Lietal., 2006, 2009] especially for small-amplitude geophysical signals with long 
wavelengths, including interseismic deformation and some anthropogenic processes [e.g., Zebker and Rosen, 
1997]. Satellite and ground-based measurements show that over land it may exhibit large spatial and temporal 
variability [e.g., Kumar et al., 201 0]. A large part of this variability may be related to the airmass transport or 
evaporation from open water bodies and evapotranspiration of vegetation and soils. In this paper, we report new 
observations pointing to a mechanism responsible for periodic oscillations in the column water vapor (CWV) 
related to mountain lee waves. 

Atmospheric gravity waves are generated by lower atmospheric sources, e.g., flow over mountains, sporadic 
diabatic heating in convective systems, and imbalance in jet streams and fronts [Kim et al., 2003]. Mountain 
lee waves are caused by an air flow over mountain ridges within a stably stratified atmosphere [Smith, 1 976; 
Durran, 1986]. Breaking waves and small-scale waves can be a source of turbulence and strong vertical air 
currents, which can be an aviation hazard [Uhlenbrock et al., 2007; Sharman et al., 2012]. Depending on 
temperature and available moisture, clouds may form in the lee of mountain ranges in rows quasi-parallel to 
the terrain and orthogonal to the direction of the flow. Satellite classification of orographic clouds began with 
the first weather satellite, the Television Infrared Operational Satellite, which was launched in 1 960 [e.g., Conover, 
1 964]. Since then, many studies reported mountain wave cloud signatures on Earth [e.g., Fritz, 1 965; Ernst, 1 976], 
and even Mars [e.g., Kahn and Gierasch, 1982], 

Atmospheric gravity waves have been widely studied using high-resolution remotely sensed temperature 
observations from space. Aqua satellite observations have been used to image and study orographic gravity 
waves in the stratosphere [e.g., Wu and Zhang, 2004; Eckermann et al., 2007; Alexander and Barnet, 2007; 
Alexander et al., 2009; Alexander and Teitelbaum, 2011]. Limb Scanning measurements from space also provide 
observations of stratospheric mountain waves [e.g., Eckermann and Preusse, 1999; Jiang et al., 2002; Alexander 
et al., 2008], Mountain-generated gravity waves from the Andes observed by the Atmospheric Infrared Sounder 
have horizontal scales of -1 0 2 km and altitudes of -20-40 km [Alexander and Barnet, 2007; Jiang et al., 2013], 
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Mountain-generated waves are also regularly observed from radiance in 6.7 (im (water vapor) channel with 
most radiance typically generated at altitudes above approximately 550 hPa [ Uhlenbrock et al., 2007; Feltzet al., 
2009] where the waves can present aviation hazards. 

The high spatial resolution (~1 km) CWV is a remote sensing product of the Moderate Resolution Imaging 
Spectroradiometer (MODIS) [Gao and Kaufman, 2003] and Medium Resolution Imaging Spectrometer [ Bennartz 
and Fischer, 2001 ; Lindstrot et at., 201 2], In this work, we are using column water vapor retrieved as part of the 
Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm developed for MODIS. MAIAC is a 
new-generation MODIS algorithm which uses a time series analysis and processing of groups of pixels for 
simultaneous retrievals of atmospheric aerosol and land surface reflectance properties [Lyapustin et at., 2011, 
201 2]. It includes cloud detection [Lyapustin et at., 2008] as well as a CWV retrieval (based on near-infrared (NIR) 
channels at 0.94 pm) required for correction of the water vapor absorption. Despite the fact that CWV is an 
operational MODIS product (MOD05) [Gao and Kaufman, 2003], MAIAC also derives CWV because the retrieval 
algorithm is fast and simple and such an approach helps avoid unnecessary data transfers and cross-product 
dependence. The MAIAC CWV retrieval approach is described in section 2, followed by a discussion of the 
mountain waves observed in CWV field in section 3. 


2. MAIAC Water Vapor Retrieval 


MAIAC retrieves column water vapor over the land surfaces and the inland or coastal waters with glint. The 
algorithm uses three MODIS NIR water vapor channels with the following band center and width (nm): 1 7 
(905, 30), 1 8 (936, 1 0), and 1 9 (940, 50). The absorption is highest in band 1 8 and decreases in bands 1 9 and 
1 7. The retrievals are based on the following approximate formula for the top of atmosphere reflectance: 

= i PAM a,M,<p), 0) 

1 - <7a0“o)c<u 


where (u 0 , /. i , <p) are cosines of solar and view zenith angles and relative azimuth, is a path reflectance, 
TfCw) is a total two-way atmospheric transmittance, c Qy is spherical albedo of atmo- 

sphere, and p A (ju 0ll u,cp) and c^tAo) are surface reflectance and albedo. The terms and T* ot are spectrally 
integrated with the spectral response function of a given channel and solar irradiance. Due to strong water 
vapor absorption and usually low aerosol scattering in the NIR, the path reflectance and the multiple scatter- 
ing of light between the surface and the atmosphere (denominator in the second term) can be omitted. If the 
surface reflectance changes little in a narrow absorption interval of 0.9-0.94 pm, then a two-channel ratio 
algorithm can be used to derive CWV: 
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For a given view geometry, the solution is found by searching the look-up table (LUT) of transmittance ratio. 
The LUT was computed with Interpolation and Profile Correction (IPC) algorithm [Lyapustin, 2003] designed 
for fast and accurate radiative transfer computations in absorption bands with arbitrary spectral resolution. 
The monochromatic water vapor absorption was computed with resolution of 0.01 crrT 1 based on HITRAN- 
2000 [Rothman et ai., 2003] using the standard atmospheric profiles. To account for different effective absorption 
and sensitivity of different band ratios, we follow [Goo and Kaufman, 2003] and compute the mean water vapor 
as follows: 

W = f,W,+f 2 W 2 . (3) 

where 1/1/ are water vapor values derived from different channel ratios and / are weighting functions related 
to the sensitivities of the two band pairs. They are computed numerically from the ratios W)/ 

7T 9 *(a 0 ,A, and Oao-A, W)/Tu(Mo,M, W) stored in the LUT. 

The developed algorithm was validated against Aerosol Robotic Network (AERONET) [Fiolben et at., 1998] CWV 
ground measurements for 156 stations globally which showed that in cloud-free conditions the retrievals are 
generally unbiased and accurate to 5-10%. These numbers agree with accuracy assessments of operational 
MOD05 product. Similarly, we found lower accuracy (20-30%) over areas with red iron-rich soils (e.g., Canberra, 
Australia), whose reflectance changes considerably in the 0.9-1 pm spectral region due to absorption of the 
iron compounds. On the other hand, MAIAC CWV is found to be lower than MOD05 by ~5-20%, providing more 
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Figure 1 . MODIS (top) Terra and (bottom) Aqua RGB and CWV image for the mid-Atlantic U.S. region on 2 June 201 1 . The 
difference in Terra-Aqua overpass time is 1 h 40 min. 

accurate result. This is confirmed by various MOD05 CWV validation studies against radiosonde, GPS, and 
Microwave Water Radiometer data showing the wet bias (higher values) of MOD05 product in the same range 
of ~5-20% [Liu et at., 201 3; Lu et at., 2011; Albert et at., 2005; Li et al., 2003; Prasad and Singh, 2009; Kumar et at., 
201 0]. At the same time, validation studies report a very good agreement of AERONET CWV with GPS and 
radiosonde measurements with r 2 ~ 0.95-0.99 [Liu et al., 201 3; Prasad and Singh, 2009; Kumar et al., 201 0], It is 
worth mentioning that AERONET CWV algorithm uses WV-band model coefficients generated by the IPC code 
[Lyapustin, 2003] for each individual Sun photometer filter function. 

Apart from possible spectroscopy-related differences in computing atmospheric transmittance that removes 
wet bias, the MAIAC approach is also simpler than MOD05 algorithm. The latter assumes a linear change of 
surface reflectance across 0.9-0.94 pm interval and uses two additional MODIS channels B2 (0.865 pm) and B5 
(1.24 pm) to predict surface reflectance in bands 17-19. Gao and Kaufman [2003] assessed the accuracy of the 
linear model as 2.4-3.9% for most soils, rocks, vegetation, and snow and 8.4% for the iron-rich soils. We tried 
the Gao and Kaufman [2003] five-band approach and found a higher noise, often by a factor of 2-3, than the 
selected three-band method, based on comparison with AERONET. The error may grow from using much 
wider spectral interval when the linear spectral model is not a good predictor for the surface reflectance in 
the 0.85-1. 24 pm region. 

At high aerosol loading and elevated profiles, the aerosol scattering increases the measured signal and the band 
ratio thus reducing retrieved CWV [see also Bennartz and Fischer, 2001 ; Lindstrot et al., 201 2]. Our AERONET-based 
analysis shows that high aerosol outbreaks explain over 90% of cases when retrieved WV was significantly lower 
than the AERONET value. These data show the need for aerosol correction on hazy days which is currently not 
implemented in either MAIAC or MOD05. 

3. Mountain Waves in MODIS CWV 

An example of mountain waves in MAIAC NIR CWV is shown in Figure 1. The Terra (top) and Aqua (bottom) 
images cover an area of 400 x 400 km 2 for 2 June (day of year 1 53) 201 1 . The MODIS RGB images (left) show 
clear conditions with low cloudiness. The average MAIAC aerosol optical depth at 0.47 pm was less than 0.2. 
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Figure 2. MODIS Aqua RGB and CVW image for 1 8 March 201 1 . The color scale is the same as in Figure 1 but with the maximal 
value of 2.5 cm. 


The right images show the retrieved CVW; the black (fill) values correspond to either dark water or detected clouds. 
The generated wave train roughly parallel to the Appalachian mountains is well visible in the lower-middle part of 
the image with a dissipation distance of -1 50 km. Weaker waves with a larger period can be seen in the northeast 
section of the Terra CWV image. The Aqua image was obtained 1 h and 40 min after the Terra overpass, so the 
wave train appears to be stationary. An apparent reduction in contrast in the Aqua CWV relative to the Terra CWV is 
a consequence of the view geometry, with an Aqua view zenith angle of -40° as opposed to the nadir view of 
Terra. In the area of maximal contrast, Terra CWV waves have an average wavelength of ~6 km with an amplitude 
of -1.0 cm, with CWV ranging from -0.5 cm in the trough to -1.5 cm in the crest. 

Figure 2 shows another example from Aqua during 1 8 March 201 1 with wavelength -7 km, in which cloud 
formation on the crests of the waves can be seen in the RGB image on the left. Finally, Figure 3 shows two 
more examples of CWV waves with very different wavelengths, in the range of 3-4 km on 3 June (top) and 
1 0-1 5 km on 1 5 October (bottom) of 201 1 . 

Generally, mountain-generated lee waves can be frequently observed east of the Appalachian range dur- 
ing winter, spring, and late fall in the MODIS NIR CWV imagery and on some summer days with low to 
moderate humidity. On such days, the air of a dry colder upper layer is alternately lifted and lowered by 



Figure 3. MODIS Aqua RGB and CWV images for (top) 3 June and (bottom) 15 October 201 1. The color scale is the same as 
in Figure 1 but with the maximal values of 2 cm (top) and 2.5 cm (bottom). 
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gravity waves, causing convergence and 
divergence of air in the warmer moist 
bottom layer, which is seen as bands in the 
CWV field. 

Atmospheric profiles of temperature and 
wind fields were used to assess the ability 
of the atmosphere to sustain trapped 
gravity waves. The vertical profile of the 
Scorer parameter (/.) calculated as N/U 
[Smith, 1976] from radiosondes and from 
the GEOS-5 analysis [Rienecker et a!., 2008] 
fields for the 2 June case (Figure 1 ) is shown 
in Figure 4. Here N is static stability and U is 
wind speed approximately normal to wave 
phase lines. The Scorer parameter [Scorer, 

1 949] arises from the dispersion relation for 
atmospheric gravity wave solutions to the 
Taylor-Goldstein equation and is used as an 
indicator of trapped lee waves. When L decreases strongly with height, trapped lee waves of certain fre- 
quencies are possible. The vertical dashed line in the figure is the horizontal wave number {k = 1 /wavelength) 
of the wave motions seen in Figure 1 , and the level at which L descends below the dashed line gives the 
trapping height of the gravity waves. Below this level ( L > k), the waves are trapped or propagate vertically, 
and above this level (L < k) the buoyancy-restoring forces cannot sustain oscillations at those high intrinsic 
frequencies and the wave motions decay. In addition, the unstable stratification in the atmospheric boundary 
layer near the surface (the region where the red line is dotted) also cannot sustain gravity waves. The trapped 
waves, therefore, based on these atmospheric conditions, can exist in a layer between approximately 1 and 
2 km above the surface, where enough water vapor exists to reveal oscillations in the CWV of MODIS. 
Although the waves decay below 1 km according to the analysis, they still penetrate some distance into the 
moister surface layer. The net result is a series of columns of alternately moister and drier air 1 -2 km deep 
near the surface that are detected by the MODIS sensor. Similar results were found for the other events in 
Figures 2 and 3, although trapping in a somewhat deeper layer below 3 km was indicated in the October case 
with the longer horizontal wavelength. 

These waves are not observed on typical summer days when CWV values of 3-5 cm are ubiquitous for the 
East Coast of the U.S. The absence of MODIS CWV observations of mountain lee waves in this season is likely 
related to deeper unstable PBL heights and atmospheric conditions that will not sustain gravity waves. 
Higher CWV would also increase the chance of forming clouds that can obscure the waves. 

4. Conclusions 

Mountain lee waves are an important phenomenon in tropospheric meteorology with practical implications 
(e.g., aviation hazards). This paper reports the first observation of lee waves in the MODIS NIR column water 
vapor east of the Appalachian mountains at high 1 km resolution. Generation and vertical propagation of 
these waves requires a stably stratified atmosphere with stability and wind conditions leading to trapping of 
shorter horizontal waves. Both radiosonde and atmospheric analysis profiles established that conditions for 
lee wave trapping below 2-3 km altitude existed on the days when the oscillations in the MODIS CWV were 
observed. Importantly, in contrast to previous studies based on sounding data which reported gravity waves 
at altitudes of the midtroposphere and stratosphere, the CWV shows perhaps the lowest near-surface waves 
that can be detected from remote sensing data over land. Knowledge of variability in CWV at these scales will 
inform regional and global models of moist processes, which rely on assumptions about the subgrid scale 
probability distribution functions for total water [e.g., Molod, 2012], Remote sensing measurements of 
greenhouse gases, such as C0 2 and CH 4 , will also benefit from improved characterization of CWV variability 
because such information allows for more accurate representation of the line broadening by water vapor and 
potential error reduction in such observations. 


IAU June l - scorer F’arameter 



L (km'’) 


Figure 4. Radiosonde (black and blue) and GEOS-5 (red) derived pro- 
files of the Scorer parameter (U2tc) at IAD on 2 June. The dotted portion 
of the GEOS-5 profile indicates the unstable planetary boundary layer 
(PBL). The dashed line indicates the horizontal wavenumber of the 
waves observed in MODIS CWV on 2 June. 
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